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af ACiLitys MATH LIBRARY 
“ABSTRACT: 
MTHSHCOSH is a function which returns the H ane a, hyperbolic cosine 


of its H floating point argument. The call is standar 
call-by-reference. 


VERSION: 1 
it Mk ede 

ed John A. Wheeler, 12-Sep-1979: Version 1 
MODIFIED BY: 
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1-006 — eet 11:24: 96 YRTHRTL SRC JMTHHCOSH.MAR; 1 . 
-TITLE MTHS$HCOSH ; 4 Shee ting Hyperbolic Cosine routine 
IDENT /1-006/ : File: MTHHCOSH.MAR EDIT: RNH1006 


1-€ 


MTHSHCOSH ; 4H 
1-006 HIS 


f 
TO 


a ting H 


Detatt 
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erbolic Cosine routine 16-SEP-1984 01: a7 08 nly Be Macro V04-00 Page 


ed Current Edit History 6=-SEP-1984 11: 


MTHRTL.SRCJ MTHHCOSH.MAR; 1 


-SBTTL HISTORY ; Detailed Current Edit History 


Edit History for Version 1 of MTHSHCOSH 


1-006 


Adapted from a PR \e, getrre 
Use MTHSH ode fy tre 4-0ct-1 
H_0.25 should be Literal rather than own storage. SBL Agg 
Don't store reserved operand before signal. SBL 7-F hy 
Changed lower Limit for Chebyshev approximation from 2**-56 to 
he Be plang second call to ~~ for » neue ett between 225, gre 

by computing COSH(x) = (Z + h Z = EXP(ix!). 
Se ae call to EXP for Reds. values between $7 ei nd 
an 
Changed att final floating point divisions by 2 to interger 
subrtracts of 1 from the exponent field. 
Extended maximum range to 16384*Ln 
cuanege’ Lo se for computing EXP tai=tn2) to reduce error. 


iliwlnetss. symbolic short Literals. RNH 15-Oct-81 


by 
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MTHSHCOSH ;_H floating Hyperbolic Cosine routine 16-SEP-1984 01: AX/VMS Macro V04-0 P 
1-006 BECLARATIONS  Goslaretive Part of Modul 6=-SEP-1984 11: 4 rk LOANS Re TATIINCOST.MAR: 1 re ad 
0090 if ~SBTTL DECLARATIONS ; Declarative Part of Module 
2000 fe} INCLUDE FILES 
Maite 6 3 
0000 77 ; EXTERNAL SYMBOLS: 
0000 78 L GBL 
44 79 EXTRN MTHSSSIGNAL 
000 80 «EXTRN MTHSK_FLOOVEMAT 
4 ee -EXTRN MTHSHEXP_R6 
0 : EQUATED SYMBOLS: 
$000 (Ba 
0000 85 
0000 86 ; 
0000 87 ; MACROS: none 
0000 88 ; 
0000 89 
0000 90 
0000 91 ; 
464 36 : PSECT DECLARATIONS: 
00000000 94 -PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
0000 95 3 program section for math routines 
0000 96 ; 
0000 97 ; OWN STORAGE: none 
0000 98 ; 
0000 99 ; 
0000 100 ; CONSTANTS: 
0000 101; 
0000 196 
0000 103 H_16383_L0G_2: 
ESEQA4SE 62DE400E 0000 104 -QUAD “XE3EQA4SE62DE400E 
6B251812 €66964DB 0008 £105 -QUAD *%*X6B251812C66964DB 3: 16383*Ln2 
0010 106 4.16384 _L0G 2: 
A39E 2F EF SESAOOE 0010 107 UAD a AS OC SER SGESS OO 
07E66730 93C7F357 0018 108 “QUAD *X07E66730935C7F 357 + 16384*Ln2 
0020 109 H_LOG_2_HI: 
asorer ee 6364000 0020 110 QUAD MAE eT ar oge sh 000 
40006730 93C7F357 43 5 4] Oe he “400067 093C7F 357 : (high 109 bits of ln2)+2*-109 
FOCB950B CODOBF9D 0030 118 “i dpe £0cB92 OBCODOBF 9D 
5F359D027 D674CD98 0038 114 - QUAD NSE 9027D674CD98 z ln2 - H_LOG_2_HI 
pp48 115 H_2_POWER_57: 
00000000 0000403A 4 116 -GUAD %X000000000000403A 
00000000 00000000 0048 £117 -QUAD 0 3: 257 
it 118 
5 119 HCOSHTAB: 
2B3EE17 acer ras Tt 120 LO F3FA8, 0 AZB SEE ITS > 2.502658279567823990162458156254711E 
ADEB4F5 13 986 005 121 $ 986, *XADEB4F5 
DO1949AF F2C935FB1 0060 1 § FB1,*XDO1949AF ; 1.611665395066462480330719767236559E 
C6034 1DF 6388C bose 1 § B8C ,*XC60346B9 
4444 D pce F 0 : 124 FBB,.“X64FF997D ; 8.896792721476872237647982378253568E 
CB38AA4 O8 F4OF 7 125 F4OF3,°X7CB38AA4 
FGF ast . 423FC 1 § 23FC + atece BAG : 4.110317621712466937445902717251860E 
$010 11 aed. 1 FAGF 1 21801 11 
9288638 68273FCC 00 128 73FCC 288638 ; 1.5619206968599445287218051911 50865E 


| 
} 
| 
| 


4.779477332387377634901255174808471E -14 
“9 
-7 
=§ 
-3 
-2 
1 
-1 
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1.147074559772972474513234083768590E -11 
2.087675698786809897912112497212423E 
2.755731922398589065255749092851911E 
2.480158730158730158730158708693922E 
1. 388888888888888888888888888904959E 
4. 166666666666666666666666666666603E 
5 .000000000000000000000000000000000E 
9.999999999999999999999999999999999E 
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MTHSHCOSH ;_H floating Hyperbolic Cosine routine 16-SEP-1984 01:35:0 AX/VMS Macro V04-00 Page 
1-006 MTHSHCOSH : W Floating point COSH funct 6-SEP-1984 91:32:88 EMTHRTL. SREIRTHHCOSH. MAR: 1 . 
} 13 -SBTTL MTHSHCOSH = H Floating point COSH function 
1 154 
1 155 ;+4 
} ! § 3; FUNCTIONAL DESCRIPTION: 
} 138 3; HCOSH = double precision H floating point function 
01 160 ; HCOSH(X) is computed as: 
0130 161; 
01 16¢ 3 If ixi-< 2*=57, HCOSH(X) = 1. 
4 165 ; IF “38° =< iXi <.0.25, HCOSH(X) = Coepyeter series 
1 164 ; If =< iki < 5 eine et Z = HEXP(iX!) and set HCOSH(X) = (2+1/2)/2 
0130 165; If 57#ln2 =< 1X! 16383(n2, then HCOSH(X) = HEXP(!X!)/2. 
0130 166 ; If 16 85s \ng =< 1X! < 16384ein2, then HCOSH(X) = HEXP(!X!=Ln2). 
61 8 19¢ 3 If 16384*ln2 =< 1X}, then overflow. 
0130 19? ; CALLING SEQUENCE: 
0130 170; 
0130 171; hcosh.wh.v = MTH$HCOSH(x.rh.r) 
0130 We : 
$138 173 ; -or- 
0130 174; 
8135 We? : CALL MTH$HCOSH(hcosh.wh.r, x.rh.r) 
0130 177; Because an H-floating result cannot be expressed in 64 bits, it is 
0130 178; returned as the first argument, with the input parameter displaced 
$138 44 ; to the second argument, in accordance with the Procedure Calling 
0130 181; 
0130 133 ; INPUT PARAMETERS: 
0130 18 
00000004 0130 184 LONG = 4 : define Longuord multiplier 
00000008 b1%6 133 x = 2 * LONG ; Contents of x is the argument 
$130 187 ; IMPLICIT INPUTS: none 
0130 188; 
$130 189 ; OUTPUT PARAMETERS: 
0130 190 ; 
00000004 0130 191 hcosh = 1 * LONG 
B15 136 ; 
0130 193 ; IMPLICIT OUTPUTS: none 
0130 194 ; 
81 8 195 ; COMPLETION CODES: none 
1 196 ; 
3 3 3 3; SIDE EFFECTS: 
6} Q 199 ; Signal: MTH$_FLOOVEMAT if 16384*Ln2 =< {Xi with reserved operand in RO/R3 
130 00 ; (copied to the signal mechanism vector CHFSL_MCH RO/R1 by IBSSIGNAL). 
1 8 01 ; Associated message is: ‘FLOATING OVERFLOW IN MATA LIBRARY’. Result is reserved 
: 4 8 3; operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1 
130 Sz ; NOTE: This procedure disables floating point underflow, enables integer 
: 0 05 ; overflow. 
130 8 
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MTHSHCOSH ;_H floating Hyperbolic Cosine routine 16-SEP-1984 01:35:0 AX/VMS Macro Vv04-00 Page 
1-006 MTHSHCOSH : WF loat ng point COSH funct etsy 01:32:88 EMTHRTL. SREIRTHHCOSH. MAR: 1 . 
1 09 
40FC 1 10 -ENTRY MTHSHCOSH, “M<IV, R2, R3, R4, RS, Rb, R7 > 
1 11 ; standard call-by-reference entry 
1 \§ ; disable DV (and FU), enable IV 
} 1 MTHSFLAG_JACKET ; flag that this is a jacket procedure in 
6D Q0000000'GF 9E 1 MOVAB G*MTHSSJACKET_HND, (FP) 
1 3; set handler address to jacket 
4 : ; handler 
0139 14 ; case of an error in routine 
6D QO0000000'GF 9E 4 9 15 MOVAB G*MTHSSJACKET_HND, (FP) ; Set handler address to jacket handler 
140 16 ; If an error, convert signal to user PC 
8128 17 ; and resigna 
50 08 BC 70FD 149 18 MOVH ax(AP), RO ; RO/R3 = {X! = @value(AP) 
50 8000 8F AA 014 19 BICW2 #*x8006, RO ; RO/R3 = 1X} 
SFFF 8F 50 B11 O14A 0 CMPW RO woXSEFF 3 compare {Xi with 0.25 
1— 18 O14 BGEQ  GEd_T0_0.25 : branch if ix! >= 0.25 
0151 $ 
0151 4; ix! < 0.25 
Bal $83 
3FC8 BF 50 B61 O151 2 CMPW RO, #*X3FC8 3 compare :Xi with 2*-57 
06 18 0156 8 BGEQ GEQ_TO_2m57 3 branch if {Xi >= 2*-57 
0158 $3 
0158 0; 
0158 31 5 IX! < 2-57 
0158 3 3 
3128 3 
04 BC 08 70FD 015 234 MOVH #1, @hcosh(AP) ; Store 1.0 as result 
04 0150 2 RET 3 return 
015E 6 
O15E * 3 
Bize 8 ; 2e-57 =< IX! < 0.25 
15E 39 ; 
015E 40 
Bize 41 GEQ_TO_2m57: 
50 50 64FD 015 $6 MULH2 RO,RO ;Get ARG*2 for POLYH. 
FEE7 CF OD 50 75FD 016 4} POLYH RO, #HCOSHLEN-1, HCOSHTAB 
016 44 3 RO/R3 = SUM(Ci*X*i) 
04 BC 50 70FD 0169 45 MOVH RO, ahcosh(AP) ; Store result in first argument 
04 OQ16E 46 RET 3; return 
oF 47 
16F re 3 
16F 49 ; 0.25 =< ix! 
16F 50 ; 
16F 2) 
16F ; GEQ_T0_0.25: 
FE8B CF 50 71FD 1¢r 5 CMPH RO, H_16383 496 2 3 compare iXi with 193850 ne 
20 14 125 34 BGTR  GTR_TRAN_16383_COG_2 =: branch if ix! > 16383*Ln2 
177 2$ 3 
177 ; 3 0.25 =< IK! =< 16383*Ln2 
17h $89 ° 
00000000 'EF 16 177 60 JSB MTHSHEXP_R6 3; RO/R3 = HEXPCiX!) 


— eee) 


> 16383*Ln2 + LOG(2) =< IX!, error 


condition value 

RO = result = reserved operand -0.0 
goes to RO/R3 so error handler 

can modify the result. 


ERROR: MOVZBL #MTHSK_FLOOVEMAT, (SP) 
ASHQ #15, #T, RO 


signal error and use real user's PC 
Restore result 
return 


ae 
CALLS #1, MTHSSSIGNAL 
mMOVO RO, @hcosh(AP) 
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MTHSHCOSH ; H floating Hyperbolic Cosine routine 16-SEP-1984 01:35:0 AX/VMS Macro V04-00 Page 
1-006 MTHSHCOSH 3 WPtosting point COSH funct aes 94:32:95 UMTHRTL. SREIMTHHCOSH. MAR: 1 . 
FEBD CF 50 71FD 017 61 CMPH RO, H_2_ POWER_57 : Compare HEXP(iX!) with 2*57, if 
9 4 66 BGTR ONE_TERM ONLY ; lar sr. only one term is needed. 
54 08 0 67FD 01 6 DIVH3 RO, #1, R4 ; R4&/R7 = HEXP(-iX!) 
50 4 60FD Q18A 64 ADDH2 =R4, RO 3 RO/R3 = HEXP(X) + HEXP(-x) 
13 65 ONE_TERM_ONLY: 
50 01 A2 Q1 ; ¢ SUBW2 #*x0001, RO 3 RO/R3 = (HEXP(X) + HEXP(-X))/2 
04 BC 50 7DFD 019 6 MOVO RO, @hcosh(AP) ; move COSH(x) to first argument 
04 8138 68 RET ; and return 
19 $2 
0197 0; 
197 71 5; 16383*Ln2 < ix 
197 ie : 
197 7 
B13% 74 GTR_THAN_16383_L0G_2: 
FE73 CF 50 71FD 0197 75 CMPH RO, H_16384_106_2 ; compare {Xi with 16384*L|n2 
1A 18 b13e 6 BGEQ ERROR ; branch to ERROR if 16384*Ln2 =<!Xx! 
Ret 78 
19F 79 ; 
019F 80 ; 16383*ln2 =< IX! < 16384*Ln2 
019F 81; 
019F Ht 
019F 8 SUBH2 H_LOG 2_HI, RO ; RO/R3=!Xi-(hi 109 bits of ln2+2*-109) 
O1AS 284 JSB MTHSHEXP_R6 > RO/R3=HEXP(!X!-H_LOG 2_HI) 
01AB MULH3) H_LOG 2_[0, RO, R4 ; R4/R7=H_LOG_2_LO#HEXP(TXi-H_LOG_2_L0) 
g18¢ ADDHS =R&, RO, ahcosh(AP) ; Store HEXP(TXT=ln2) in 1st argument 
018 RET 3 return 
01B9 
01B9 
01B9 
0189 
01B9 
01B9 
01BD 
01C1 
01C1 
01C1 
01C¢3 
O1CA 
O1CF 
01D0 
01D0 
01D0 
01D0 +3 
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; H floating Hyperbolic Cosine routine 
Symbol table 


ERROR 6° 1B9 R 1 
GEQ_T0_0.25 OOO16F R 1 
GEQ_TO_2m57 0001 5 R 1 
GTR-THAN_16383_L0G_2 0001 R 1 
HCOSH z 3 4 
HCOSHLEN = 0 f 

OSHTAB B86 R 01 
W1038)-L06_ 000000 R 01 
H_16384_ LO it 10 R 01 
H_2_ POWER 57 0000040 R 4 
WLOG ¢HT 000 $8 R 1 
H_LOG_2_L0 $38 & 01 
LONG = 00 04 
MTHSSJACKET_HND eerkeree x 4 
MTHSSSIGNAL reekeree x 0 
MTHSHCOSH 00000130 RG 01 
MTHSHEXP_R6 eeekeene x 00 
MTHSK_FLOOVEMAT tererere xX 00 
ONE_TERM_ONLY 0000018E R 01 

= 00000008 
pmenmen awe wrnnraoceae + 
! Psect synopsis ! 

PSECT name Allocation PSECT No. Attributes 
es, 00000000 <¢ 0.) 00 ¢ OO.) NOPIC USR CON ABS) LCL 
_MTHSCODE 000001D0 (¢ 464.) O01 ¢ 1.) PIC USR CON REL LCL 


ee 


Phase Page faults CPU Time Elapsed Time 
Initialization 32 Bb $44 14 00:00:01.11 
Command processing 126 00:00:00.68 99:88:85 -06 
Pass 1 91 00:00:01.06 0:00:04.74 
Symbol table sort 0 B88: 0098 89 288 : 20-00 
Pass 2 66 0:00:00. 00: Be Ae 
Symbol table output 4 00:00: 8-05 00:00:0 “98 
poses synopsis cutout 3 Bo 282: a: 09 :00:00.06 
ross-reference ou :00:00. :00:00. 
Assembler run fetal 324 $080; 3°88 8:96: 90-98 


The working set Limit was 900 pages. 
shez bytes +) pages) of virtual memory were used to buffer the intermediate cod 
ere were 
364 source Lines were read in Pass 1, producing object records in Pass 
1 r xge of virtual memory was used to define 1 macro. 


ESEp=1986 11:30:69 


e. 
pages of symbol table space ah locates to hold 20 non-local ye 0 local symbols. 


AX/VMS Macro V04-00 Page 
MTHRTL.SRCJMTHHCOSH.MAR; 1 


-— NOEXE NORD NOWRT NOVEC BYTE 


EXE RD NOWRT NOVEC LONG 


pee aa : a 
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VAX=-11 Macro Run Statistics — n “SEP=1 1382 9: 32 08 MTHRTL.SRC JMTHHCOSH.MAR; 1 ° (4) 


$ eer ere ereceann een r wee cece eeeece + 


: Macro Library Statistics : 


Macro Library name Macros _defined 
“S255SDUAZ8:(SYSLIBISTARLET.MLB;20 0 itt tstCtCtCS~S~S 5 
O GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHCOSH/OBJ=0B/$:MTHHCOSH MSRC$:MTHJACKET/UPDATE=(ENH$:MTHJACKET) #MSRC 
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